Monopole confinement by flux tube 
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Abstract 

\Q \ We revisit Nambu's model of quark confinement by a tube of magnetic flux, with two additional 

features. The quarks are taken to be magnetic monopoles external to the tube, which seal the 
ends, and are also taken to be fermions. This ensures that the model is inconsistent unless there 
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are at least two species of fermions being confined. 
PACS numbers: 11.27.+d, 14.80.Hv 
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Many years ago, Nambu suggested a model for the confinement of quarks [1], in which 
quarks are the endpoints of an open Dirac string [2] . If the gauge field is massive, this string 
becomes real [3], unlike the usual Dirac string, which is a gauge artifact. In this paper 
we make a small modification to Nambu's model - instead of treating the endpoints as 
quarks, we consider quarks as external fermionic monopoles, which are then used to seal the 
endpoints of the open string. This ensures that flux does not escape or enter the endpoints, 
and a delta function(al) ensures that quark currents always stick to the ends of strings. The 
string itself is for us an infinitesimally thin tube of magnetic flux, created in an Abelian 
Higgs model. 

In the presence of strings, we first dualize the scalar field to find the strings interacting 
via an antisymmetric tensor potential [4, 5], while the Abelian gauge field is dualized to a 
'magnetic' photon [6, 7]. Next we introduce fermionic magnetic monopoles into the theory 
and minimally couple these to the magnetic photon. The monopole current has to be 
axial, which produces an anomaly. We cancel the anomaly by postulating additional species 
of fermionic monopoles. Then we dualize the resulting theory again, to find a theory of 
magnetic flux tubes interacting with a massive Abelian vector gauge field. The tubes are 
sealed at the ends by fermions, thus providing a toy model for quark confinement. 

We start with the generating functional for the Abelian Higgs model in 3 + 1 dimensions, 
coupled to an Abelian gauge field A e . The partition function is given by 



S = J d 4 x (-^F^F"" + \D^\ 2 + A(|$| 2 - v 2 f^j , 



(1) 



where = <9 M + ieA^ , and F^ v = d^A e v — d^A^ is the Maxwell field strength. 

In order to consider the situation where magnetic flux tubes are present, we change 
variables from $ ,$* to the radial Higgs field p and the angular field 9 , defined by $ = 
-^pexp(i^) . Then the measure becomes, in these variables, 



/ 



mm* ■■■= 9p 2 m ■■■ . (2) 



where the dots represent the measure for any other fields and the integrand. Here we will 
consider the theory in the London limit (A — > oo with the radial part of the field fixed), i.e. 
we will ignore the p-dependent part of the measure, and set p — v (constant) in the action. 

In the presence of flux tubes we can decompose the angular field into a regular and a 
singular part [8], 9 = 9 r + 9 s , where 9 s corresponds to a given magnetic flux tube, and 9 r 
describes single valued fluctuations around this configuration. The singular part of the phase 
of the Higgs field is related to the world sheet E of the magnetic Abrikosov-Nielsen-Olesen 
string according to the equation, 

e^ px d p d x 9 s = IT, (3) 
JT = 2nn j da^(x(0) S\x - x(0) , (4) 

where £ = (£ x ,£ 2 ) are the coordinates on the world-sheet of the flux-tube, and da pu (x(^)) = 
e ab d a x tl d} ) x 1 ' . In the above equation 2n is the vorticity quantum in the units we are using 
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and n is the winding number [9]. The partition function now has the form 

1 „,2 



A* WW exp 



d A x[--F^F^ + -(d,9 + eAl 



e\2 



(5) 



We can linearize the term — (d p 6 + eA p ) by introducing an auxiliary field C p to get 
J W exp y d 4 x^(d p 6 r + d p 6 s + eAlf 



= J WW M exp 
= J @C„5[d^} exp 



-i J d'x^ 2 Cl + C^d,9 s + eA;) 



(6) 



We can resolve the constraint <9 M C M = by introducing an antisymmetric tensor field B pv 
and writing C p in the form C M = \e pvpX d v B p \. Integrating over the field C p , we get 



Z = J ^A e ^x p ^)&B pu exp i J d 
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■ l —H^ p H^o - hl, v B^ - e -e^ x Ald v B i 



(7) 



Here we have written H pvp = d p B up + d u B pp + d p B pv . We have also replaced the integration 
over *3)d s by an integration over *3)x p {£) which represents a sum over all configurations of 
the worldsheet of the flux tube. Here x p {£) parametrizes the surface on which the field 9 
is singular. The Jacobian for this change of variables gives the action for the string on the 
background spacetime [8] . The string has a dynamics given by the Nambu-Goto action, plus 
higher order operators [10], which can be obtained from the Jacobian. We will not write the 
Jacobian explicitly in what follows, but of course it is necessary to include it if we want to 
study the dynamics of the flux tube. 

Let us now integrate over the field A e . To do this we have to linearize F^ p F e ^ u by 
introducing another auxiliary field x^v , 



J W^exp i J d 4 : 

S^W^exp 



.I F e pep, _ - e ^P^A e u d u B 



d x 



\x^ v + \^ P \,ud p Al - 6 -e^B pu d p Al 
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e^d^XpX ~ eB pX ) exp[* / ^{--x^xT}] 



We can integrate over Xnv by solving the 5-functional as 

X P u = eB pu + d u A 



'If >- v 



(8) 



(9) 



Here A™ is the 'magnetic photon', and what we have achieved is dualization of the vector 
potential A^. The result of the integration is then inserted into Eq. (7) to give 



Z = J ®A™®x,(0®B„ u e W i J L^eB^ + d^f + ^H^H^-^B^ 



(10) 
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The equation of motion for the field can be calculated from this to be 



TP ^ 

d x H x » v = G^-rrpYT, (11) 

e 

where G^ u = eB^ u + d^A™ — d v A™ , and m = ev . Absence of magnetic current gives 
= 0, so from Eq. (11) we get 

d v YT = . (12) 

In particular, this equation means that in the absence of magnetic monopoles, the vorticity 
tensor current E^j, is conserved. In other words, due to the conservation of magnetic flux 
all the flux tubes in this case are closed or infinite. The calculations so far have been done 
previously by others. What we will do now has not appeared in the literature before this. 

We now attach magnetic monopoles at the ends of a finite flux tube. We will take the 
monopoles to be massless fermions and minimally couple the monopole current to the mag- 
netic or dual photon. The monopole current must be axial for compatibility with Maxwell's 
equations. After coupling, we will dualize the theory a second time to get back to vector 
gauge fields, now coupled to flux tubes. 

However, a theory containing axial fermionic currents is anomalous and if we try to 
dualize the theory, the presence of the anomaly gives inconsistent results. We can cancel the 
anomaly by introducing another species of fermionic monopoles with axial charge opposite 
to the previous one. Let us denote the two species by q and q' , with monopole charges +g 
and — g , respectively. So the monopole current becomes 

& = gqi^q - gq'wW • (13) 

The partition function of Eq. (10) is modified to include the fermionic monopoles, mini- 
mally coupled to the 'magnetic photon' A™ , so the Lagrangian reads 

2 = -\(eB»„ + d^f + -±- 2 H, vp H^ - l^B^ + tq$q + tq'$q' - A™f m . (14) 
The conservation condition of Eq. (12) is now modified to 

-d v ^=f m . (15) 

This equation is a consequence of gauge invariance. To see this, we take a transformation 

— > B^ v + (9 M A^ — d v K^ , 
A™ - A? - - g A, . (16) 

The second term of the Lagrangian of Eq. (14) is invariant under the above transformation, 
while the first term is made invariant by setting eg = k . This is related to the Dirac 
quantization condition as we shall see shortly. 

Since the flux due to the monopoles is fully confined in the tube, the flux in the tube 
is 4ng. Using Eq. (4), Eq. (15) and the condition eg = k , we see that the "vorticity flux" 
must be equal to flux of the monopoles. So we can write 

k 

2nn = vorticity flux = —(monopole flux) = Ank (17) 
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So k = — , and we must have the Dirac quantization condition eg = — , from our assump- 

tion that the monopole flux is fully confined in the tube. We can now write the partition 
function as 



@A™@x il {£ > )@B lxv @q@q@q'@q" expi / d 



x 



1 



m\2 



- (eB^ + d^) 



+ j^H, up H^ - \^ V B^ - ZTd^ + eh A + iq(j)q + iq'$q' - A™f m 



(18) 



Since (16) is only a change of variables, Z cannot depend on A^. Thus A M can be integrated 
out with no effect other than the introduction of an irrelevant constant factor in Z, which 
we ignore. After integrating over A M , we get 



Z = I ■ ■ ■ 5 



-d^+f m exp* / d*x 



\{eB pu + d [p A^ + -^H pup H^ 



^B^ + iqfa + iqifai-A™^ 



, (19) 



where the dots represent the measures for the other fields and x^. One can see from the 
5-functional that the vorticity current tensor is not conserved, but is cancelled by the current 
of the added fermions. So the strings are open strings with fermions stuck at the ends. Now 
we dualize the theory a second time and get back to a vector gauge field. 

Introducing an auxiliary field x^v t° linearize the first term of the Lagrangian, we get 



exp i Id x 



~XiivX 



+ ^ pX x,A^x 



2 



+^ x X »vB pX + A-H^H^ - \^B^ + iqfa + zq'$q' - A™f m \ . (20) 
We can now integrate out Af^, and the result is 



@X 



i 



vXp\ J 



exp i id x 



\X\ivX 



\1V 



+ ^ upx X ^B pX -, 



l —H pup H^ - + iqfa + iq'$q'] . (21) 



1 1 

Both the 5-functionals must be satisfied, which requires -d v YT v — -t pupx d P Xp\ = . This 

can be solved by introducing a gauge field A p , which allows the integration over x^v ■ Then 
the partition function becomes 



Z = J ^x^B^A^ ■ ■ ■ 6 [i^S^ + #J exp i J d 4 x [ - -F pv F^ 
+ -^H pup H^ + ±e^ X B^d p A x + iq$q + tq'$q'] . 



(22) 



Here F^ u = d^A^ — d u A p — —e Pil/(T \L uX , the dots represent the fermion measure and we 
continue to suppress the action for the flux tube itself, as in Eq. (7). 
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The theory is now in the form we originally intended, and contains thin tubes of flux. The 
new feature is that the ends of the flux tube are sealed by fermions, so that no flux escapes, 
all flux is confined. We should not think of this as any more than a toy model of confine- 
ment, because the underlying theory is only the Abelian Higgs model and not quantum 
chromodynamics. Even then, some features are interesting enough to be highlighted. 

There is a simple argument to calculate the length of the string. The flux confined inside 
the string is Ang, a constant. The radius of the string core is of the order of 1/vy/X. From 
this we can calculate the energy per unit length of the tube to be \x ~ g 2 v 2 \ , also a constant. 
Such a string, of finite length, would collapse in order to minimize the energy unless it was 
stabilized by its angular momentum. For a rotating string of length / , energy per unit length 
/i, angular momentum J , the energy function is E — + J 2 /2{il 3 . This has a minimum for 
the length L ~ yj J j \i . We see that for the stable flux tube with magnetic monopoles at the 
ends, 

— = constant , (23) 
E 2 

similar to the well-known Regge trajectory for mesons. 

The gauge field is massive, with mass m = v/g [11, 12]. It does not couple directly 
to the fermionic monopoles at the ends. Those fermions are coupled only through the 
^-functional in Eq. (22), which guarantees that the monopoles must seal the ends of the 
string. However, any other gauge field, Abelian or not, axial or not, may be coupled to 
these fermions with charge assignments independent of their charges under A™ , which has 
been integrated out of the theory. In particular, if we suggestively rename q and q' to u 
and d, the allowed configurations are ud,ud, and uu±dd , which can couple to electroweak 
gauge fields. Note also that we could have introduced three species of fermions (with charges 
1, 1, —2, for example) in Eq. (13), for the purpose of anomaly cancellation, and we would 
again get flux tubes with ends sealed by fermions. But a single species of fermions would 
not produce such configurations. 

In order to make the model more realistic, one would need to check if flux is truly confined 
in the tube or if it escapes when the tube has a finite thickness. A similar picture starting 
with an axial gauge field and ending with a tube of 'electric' rather than magnetic flux 
will be interesting as well. Further, the freedom to have other global symmetries in the 
theory allows in principle that the U(l) producing the string here may be embedded in 
an SU(N) gi b a i x U(l)i oca i symmetry, as in [13]. We have not investigated these questions, 
and also not the million-dollar question, which is whether quantum chromodynamics has a 
low-energy effective sector that behaves like an Abelian Higgs model. 
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